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Mesoscopic wave turbulence. 
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We report results of sumulation of wave turbulence. Both inverse and direct cascades are observed. The 
definition of "mesoscopic turbulence" is given. This is a regime when the number of modes in a system involved 
in turbulence is high enough to qualitatively simulate most of the processes but significantly smaller then the 
threshold, which gives us quantitative agreement with the statistical description, such as kinetic equation. 
Such a regime takes place in numerical simulation, in essentially finite systems, etc. 
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The theory of wave turbulence is developed for in- 
finitely large system. In weakly nonlinear dispersive me- 
dia, the turbulence is described by a kinetic equation for 
■squared wave amplitudes (weak turbulence). However, 
all real systems are finite. Computer simulation of wave 
'turbulence can also be perfomed only in finite system 
(typically in a box with periodic boundary conditions) . 
It is important to know how strong discreteness of a 
isystem impacts the physical picture of wave turbulence. 

Let a turbulence be realized in a Q-dimensional cube 
■with side L. Then, wave vectors form a cubic lattice 
'with the lattice constant Ak — 27t/L. Suppose that 
four-wave resonant conditions are dominating. Exact 
resonances satisfy the equations 



k ki - k2 - kg = 0, 



(1) 



A = Lu{k) + Lu{ki) - Lu{k2) - uj{k3) = 0. (2) 

In infinite medium, Eqs. ^ and ||2Jl define hyper- 
surface of dimension SQ — 1 in 4(5-dimensional space 
'k, ki , k2 , ka . In a finite system, and (|2Jl are Dio- 
.phantine equations which might have or have no exact 
solutions. The Diophantine equation for four-wave res- 
onant processes are not studied yet. For three- wave res- 
onant processes, they are studied for Rossby waves on 
the /3-plane 0. 

However, not only exact resonances are important. 
Individual harmonics in the wave ensemble fluctuate 
with inverse time Fk, dependent on their wavenumbers. 
Suppose that all Fk^ for waves, composing a resonant 
quartet, are of the same order of magnitude Fk^ ^ F. 
Then resonant equation Q has to be satisfied up to 
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accuracy A ^ F, and the resonant surface is blurred 
into the layer of thickness 6k /k ~ Tk/oJk- This thick- 
ness should be compared with the lattice constant Afc. 
Three different cases are possible 

1. 5k ^ Afc. In this case the resonant layer is thick 
enough to hold many approximate resonant quar- 
tets on a unit of resonant surface square. These 
resonances are dense, and the theory is close to 
the classical weak turbulent theory in infinite me- 
dia. The weak turbulent theory offers recipes for 
calculation of F^. The weak-turbulent F^ are the 
smallest among all given by theoretical models. 
To be sure that the case is realized, one has to use 
weak-turbulent formulae for F^. 

2. 5k < Ak. This is the opposite case. Reso- 
nances are rarefied, and the system consists of a 
discrete set of weakly interacting oscillators. A 
typical regime in this situation is the " frozen tur- 
bulence" 121, Ell which is actually a system 
of KAM tori, accomplished with a weak Arnold's 
diffusion. 

3. The intermediate case 5k ~ Afc can be called 
" mesoscopic turbulence" . The density of approx- 
imate resonances is high enough to provide the 
energy transport along the spectrum, but low 
enough to guarantee " equal rights" for all harmon- 
ics, which is a necessary condition for applicability 
of the weak turbulent theory. 

In this article we report results of our numerical ex- 
periments on modeling of turbulence of gravity waves 
on the surface of deep ideal incompressible fluid. The 
motivation for this work was numerical justification of 
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Hasselmann kinetic equation. The result is discovery of 
the mesoscopic turbulence. The fluid motion is poten- 
tial and described by shape of surface ri{r, t) and veloc- 
ity potential ■(/'(r, i), evaluated on the surface. These 
variables satisfy the canonical equations 



dt 



SH 



'dt 



SH 
Sr/ ' 



(3) 



Hamiltonian H is presented by the first three terms in 
expansion on powers of nonlinearity V77 



H = Ho + Hi + H2 + 

grf + ) dxdy, 



Ho 



Hi 



H, 



V 



dxdy, 



(4) 



ri{kip) k{ri(k^)) + t^V^tA da;d?/ 



Thereafter, we put gravity acceleration equal to g = 1. 
Here, fc is a linear integral operator — V— V^^ , 
such that, in fc-space, it corresponds to multiplica- 
tion of Fourier harmonics ('0k — ^ J dxdy) by 
'^k'^ + ky. For gravity waves, this reduced Hamiltonian 
describes four-wave interaction. Then, dynamical equa- 
tions (PI acquire the form 

T] = kijj — (V (rjV'tjj)) ~ k[r]k^j] + 

+k{'nk[vh]) + ^"^^[v^h] + iM^y^V^v^], 



(VV)2 - (k^f 



-9V 



-[k'ip]le[r]k'il!] — [r]kip]V^^p. 



(5) 

Let us introduce the canonical variables Ok as shown 
below 



flk 



2fc'^'' + *V^^''' 



(6) 



where ujk — Vgk. In these so called normal variables 
equations 10) take the form 



dak 
dt 



. SH 



(7) 



The physical meaning of these variables is quite clear: 
|akp is an action spectral density, or jokpA/c^ is a num- 
ber of particles with the particular wave number k. 

We solved equations © numerically in a box 27r x 2n 
using a spectral code on rectangular grid with dou- 
ble periodic boundary conditions. The implicit energy- 
preserving scheme, similar to used in 0, [7], HI, was 
implemented. We studied evolution of freely propagat- 
ing waves (swell) in the absence of wind in the spirit 
of paper (^j. Different grids (512 x 512, 256 x 1024, 



256 X 2048) with different initial data were tried. In 
all the cases, we observed mesoscopic wave turbulence. 
The most spectacular results are achieved on the grid 
256 X 2048. 

As initial conditions, we used a Gauss-shaped distri- 
bution on a long axis of the wavenumbers plane 



|ak| = Ai exp 




,|k-ko| < 2A, 



|ak| = 10-^Mk-ko|>2A, 
Ai^hy. 10"^ Di = 30, ko = (0; 150). 

(8) 

The initial phases of all the harmonics were random. 
The average steepness /i =< |V7y| >~ 0.115. To stabi- 
lize the computations in the high-frequency region |10 |. 
we introduced artificial damping, mimicking viscosity 
at small scales, and an artificial smoothing term to the 
equation for the surface evolution 

90k 90k 



dt 



Ik 



(9) 



drjk 
dt 

0, fc < fcrf, 
^ -j(k~ kdf,k> fcd, 
fcd = 512, 7 = 2 X 10^ T = 3.1 X 10""*. 

With the time step r, this calculations took about two 
months on AMD Athlon 64 3500-1- computer. During 
this time, we reached 1500 periods of the wave in the 
initial spectral maximum. 

The process of waves evolution can be separated in 
two steps. On the first stage (about fifty initial wave pe- 
riods), we observe fast loss of energy and wave action. 
This effect can be explained by formation of " slave" har- 
monics taking their part of motion constants. Initially 
smooth spectrum becomes very rough. The spectral 
maximum demonstrates fast downshift. 

In the second stage, the downshift continues but all 
processes slow down. Plots of energy, wave action, mean 
frequency, and mean steepness are presented in Figs. 
I1I4I One can see clear tendency to downshift of 
the spectral maximum corresponding to inverse cascade; 
however, this process is more slow than predicted by the 
weak turbulence theory. Self-similar downshift in this 
theory gives JT], [T^ 

In our experiments 

where a decreases with time from 1/16 to 1/20. Evo- 

27r 

lution of angle averaged spectra Nk = / jakpfcdfcdz? is 
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FiglU Total energy of the system. 
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FigE] Mean steepness of fluid surface. 
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FiglU Total action of the system. 
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FigEl Averaged with angle spectra. Downshift of spec- 
tral maximum is clearly observable. 
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FigOl Frequency of the spectral maximum. 



presented on Fig. [S] Their tails (Fig. O are Zakharov- 
Filonenko weak-turbulent Kolmogorov spectra cor- 
responding to direct cascade 
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This result is robust; it was observed in similar calcula- 
tions 0, El, [g. 





t=648, <|Vri|>.0.61x10"' 

k" 































(10) 



Fig'iniTails of angle-averaged spectrum in double loga- 
rithmic scale. T — 648 = 1263ro. Power-like tail and 
front slope are close to predicted by the weak turbulent 
theory. 
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Two dimensional spectra in the initial and in the last 
moments of calculations are presented in Fig. [7| |S| One 
can see formation of small intensity "jets" posed on the 
Phillips resonant curve ^1] 



2tj(ko) = t^(ko + k) + ujiko - k). 



(11) 
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Fig|7| Level lines of logarithm of initial spectra distri- 
bution. T = 0. 




FiglHl Level lines of logarithm of spectra distribution 
at r = 648 = 1263ro. 



The spectra are very rough and sharp. The slice of 
spectra along the hue (0; ky) in the end of the computa- 
tions is presented on Fig. |^ Evolution of squared wave 
amplitudes for a cluster of neighbouring harmonics is 
presented in Fig. ^1 

Results presented in Fig. ^| show that what we 
modeled is mesoscopic turbulence. Indeed, character- 
istic time of amplitude evolution on a figure is a hun- 
dred or more their periods; thus T/lu^ is comparable 
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FigEl Slice of spectrum on axis (0; ky) at T — 648 
1263ro. 
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Figlini Evolut ion of some cluster of harmonics and a 
distant large harmonic. 



with Ak/k. On the same figure we can see the most 
remarkable features of such turbulence. 

The weak turbulence in the first approximation 
obeys the Gaussian statistics. The neighbouring har- 
monics are uncorrelated and statistically independent 
((afeafe+i) = 0). However, their averaged characteristics 
are close to each other. This is a "democratic society". 
On the contrary mesoscopic turbulence is an " oligarchic 
society". The Phillips curve I|1HI has a genus 2. Af- 
ter Faltings' proof ^21 of Mordell's hypothesis ^Hl we 
know that the number of solutions of the Diophantine 
equation 

A ^ 2(n2 + to2)1/4 _ ^ _^ (^j^ ^ 
-[(n-x)^ + (m-y)2]i/4 = q 

(12) 

is at most finite and most probably, except for a few 
trivial solutions, equals to zero. The same statement 
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is very plausible for more general resonances. Approxi- 
mate integer solutions in the case 

|A|<e 

do exist, but their number fast tends to zero at e 0. 
Classification of these solutions is a hard problem of the 
number theory. These solutions compose the "elite soci- 
ety" of the harmonics, which play the most active role in 
the mesoscopic turbulence. Almost all the inverse cas- 
cade of wave action is realized within members of this 
"privileged club" . The distribution of the harmonics ex- 
ceeding the reference level |afep — 10"^^ at the moment 
t = I2OOT0 is presented in Fig. ^] The number of such 
harmonics is not more than 600, while the total number 
of harmonics involved into the turbulence is of the order 
of 10*. 
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Fig llll Harmonics with square modulus exceeding the 
level 10"" at T = 648 = 1263To. 

Note that a situation with direct cascade is differ- 
ent. As far as the coupling coefhcient for gravity waves 
growth as fast as with the wave number, for short 
waves Tk/iOk easily exceeds Afc/fc, and the conditions 
of applicability of the weak turbulent theory for short 
waves are satisfied. 

Note also that the mesoscopic turbulence is not a 
numerical artefact. Simple estimations show that, for 
gravity waves, it is realized in some conditions in basins 
of a moderate size, like small lakes as well as in experi- 
mental wave tanks. It is also common for long internal 
waves in the ocean and for inertial gravity waves in at- 
mosphere, for plasma waves in tokamaks, etc. 
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